* Plot IRF
local labely : variable label $depvar
local yname = subinstr("`labely'"," ","",.)

local plotmin = 0
local plotmax = $hmax - 1

cap gen QuartersPlot = Quarters - 1 

local j = ${j}

//forv j = 0/1 {

	if ${cfbands} == 1 {
		local lc1 = "${lc_1cfband}"
		local lp1 = "${lp_1cfband}"

		twoway ///
			(connected b`j' QuartersPlot, lcolor(blue) lpattern(solid) lwidth(medthick) msymbol(oh) mcolor(blue) ) ///
			(line u${alpha1}`j' QuartersPlot, lcolor("`lc1'") lpattern("`lp1'") lwidth(thin)) ///
			(line d${alpha1}`j' QuartersPlot, lcolor("`lc1'") lpattern("`lp1'") lwidth(thin)) ///
			(line Zero QuartersPlot, lcolor(gs13)), legend(off) aspect(0.8, place(center)) ///
			title("{bf:Response of `labely' to a 1% of GDP tax cut}", color(black) size(medium) ) ///
			xlabel(`plotmin'(1)`plotmax' ) ytick( `yrange' , angle(0))   ///
			ytitle("Percent deviation", size(small)) xtitle("Quarters", size(small)) ///
			plotregion(margin(zero) color(white)) ///
			graphregion(color(white) margin(small)) name(IRF_`yname'_`j', replace) 

	}
	
	else if ${cfbands} == 2 {
		local lc21 = "${lc_2cfband_alpha1}"
		local lp21 = "${lp_2cfband_alpha1}"
		local lc22 = "${lc_2cfband_alpha2}"
		local lp22 = "${lp_2cfband_alpha2}"

		twoway ///
			(connected b`j' QuartersPlot, lcolor(blue) lpattern(solid) lwidth(medthick) msymbol(oh) mcolor(blue) ) ///
			(line u${alpha1}`j' QuartersPlot, lcolor("`lc21'") lpattern("`lp21'") lwidth(thin)) ///
			(line d${alpha1}`j' QuartersPlot, lcolor("`lc21'") lpattern("`lp21'") lwidth(thin)) ///
			(line u${alpha2}`j' QuartersPlot, lcolor("`lc22'") lpattern("`lp22'") lwidth(thin)) ///
			(line d${alpha2}`j' QuartersPlot, lcolor("`lc22'") lpattern("`lp22'") lwidth(thin)) ///
			(line Zero Quarters, lcolor(gs13)), legend(off) aspect(0.8, place(center)) ///
			title("{bf:Response of `labely' to a 1% of GDP tax cut}", color(black) size(medium) ) ///
			xlabel(`plotmin'(1)`plotmax' ) ytick( `yrange' , angle(0))   ///
			ytitle("Percent deviation", size(small)) xtitle("Quarters", size(small)) ///
			plotregion(margin(zero) color(white)) ///
			graphregion(color(white) margin(small)) name(IRF_`yname'_`j', replace) 
	}

//}
